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Heavy Light Weak Matrix Elements* 
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I review the status of lattice calculations of heavy-light weak matrix elements, concentrating on semileptonic B decays 
to light mesons, B—*K*'y, the B meson decay constant, /s, and the mixing parameter Bb- 



1. INTRODUCTION 

In this review I will focus on calculations of se- 
lected matrix elements involving 6 quarks. The 
semileptonic decays B^^tt^I'vi and B^^p^l~vi 
depend on the CKM element Vub which fixes one side 
of the unitarity triangle. Lattice calculations open 
the possibility for a model independent extraction 
of this quantity. The rare radiative decay B—fK*j 
is important for the determination of Vts and as a 
window on new physics. Calculations here are ham- 
pered by the necessity to impose models for form 
factors. For the B meson decay constant, the sys- 
tematic uncertainty arising from continuum extrap- 
olation is being reduced and quenching effects are be- 
ing addressed. Combining with results for the mixing 
parameter, the phenonienologically relevant quantity 
f^Bs will become better determined. 

2. SEMILEPTONIC AND RADIATIVE 
HEAVY-TO-LIGHT DECAYS 

Lattice calculations of form factors are crucial for 
these decays: heavy quark symmetry relates differ- 
ent matrix elements and hence form factors, but 
the overall normalisation at the zero recoil point 
(jj = vv' — 1 is not fixed by the symmetry as it 
is for a heavy-to-heavy transition. Here v and v' are 
the four-velocities of the b quark and the light quark 
it decays into, respectively. 

In reviewing results for these decays, I would like 
to emphasize how lattice calculations may be used 
for model independent extractions of CKM angles. 

The diagram in Fig. |^ illustrates the problem in 
obtaining form factors for B decays over the full 
range of = NP -\- — 2Mmuj, when a heavy 
meson of mass M decays to a light meson of mass 
m, transferring four-momentum q to the leptons or 
photon. Calculations using the conventional lattice 
Dirac equation (hereafter referred to as the "con- 
ventional" approach) are currently performed with 
heavy quark masses around the charm mass. The ini- 
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Table 1 

Leading M dependence of form factors for heavy-to- 
light decays in the helicity basis. The dependence 
follows from heavy quark symmetry applied at fixed 
velocity transfer lu = v ■ v' . Note that only three of 
the four Ai form factors for B'^^p^l~i'i are indepen- 
dent. 
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tial heavy meson is given or 1 lattice units of three- 
momentum, while the light final meson can generally 
be given up to two lattice units of spatial momen- 
tum, allowing g2 to be varied from g^j^^ — {M—m)'^ 
(where uj = 1) down to < at the D scale. 
Heavy quark symmetry shows that the the form fac- 
tors scale simply with the heavy mass M for fixed oj. 
Table |^ shows the leading M dependences at fixed 
LU for the form factors in the helicity basis where 
the quantum numbers are fixed for the i-channel. 
These are multiplied by power series in 1 /M to build 
up the full M dependence. 

As Fig. shows, fixed-w scaling sweeps all the mea- 
sured points to a region near g^^x f^^' ^ decays. The 
problem is then to extrapolate to cover the range 
back down to q^ ~ 0. This is particularly acute for 
the radiative decay B^K*^ where the two-body fi- 
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Figure 1. range for heavy-to- light decays as a 
function of the decaying heavy meson mass. Lines of 
constant w are shown. The hght final state mass is 
taken to be 850 MeV, typical of lattice pseudoscalar 
or vector meson masses before chiral extrapolation. 



nal state with a real photon implies that only the 
form factors at = contribute to the decay. Even 
when the heavy quark is treated using a static, non- 
relativistic or other modified action, the restriction 
on the available three-momenta on current lattices 
ensures that form factor values are obtained only 
near g^^^. 

Some assistance is provided by ensuring that any 
model dependences employed satisfy the require- 
ments of heavy quark symmetry together with any 
known constraints relating form factors at q^ = 0. 
For example, /+ and in B ^ tt decays are related 
at q^ — 0, and consistency is achieved by fitting /+ to 
a dipole [pole] form and /° to a pole [constant] form. 
Heavy quark symmetry and light flavour SU{i) sym- 
metry relate the form factors for B^^p'^l~vi and 
B~^K*j, and models can further relate these form 
factors to those for B'^^n^l^vi. An overall fit to 
all the form factors might then be used. However, 
it is clearly desirable to avoid models entirely: I will 
concentrate on model independent results below. 

Before discussing lattice results, I will summarise 
the experimental situation. The decays 5°— 
and B'^ ~>Tr^ Di are measured using full reconstruc- 
tion by the CLEO collaboration, who find [^: 

BiB'^^iT+l-ui) = 1.8(4)(3)(2) X 10-^ (1) 

BiB'^^p+l-Di) = 2.5(4)(?)(5) X 10-^ (2) 

The experimental efficiency is model dependent, so 
the above results are a combined value using different 
models. When converted to a result for \Vub\, the 
model dependence dominates the errors: 

|Kf,| = 3.3(2)(3)(7) X 10-3, 



to be compared with |Kb| = (3.1 ± 0.8) x 10"^ ob- 
tained from the lepton energy spectrum endpoint 
analysis in inclusive decays, which is also model de- 
pendent. The possibility of using lattice results to 
remove the model dependence from the determina- 
tion of \Vub\ is exciting motivation to pursue such 
studies. 

The experimental result for B—fK*j also comes 
from CLEO, who quote [|j, 

B{B^K*"f) = 4.2(8)(6) x 10"^ (3) 

Combining this with their measurement for the in- 
clusive decay ||], 



Bib -> S7) = 2.32(57)(35) x 10" 



they find 



Rk' 



T{B^K*-i) 



0.181(68). 



(4) 



(5) 



2.1. Semileptonic B ^ n 

Lattice results for the B^^n^l^vi form factors, 
/+ and /o, are available from ELC [& APE (l|Jl|], 
UKQCD and Wuppertal-HLRZ]p|P] , together 
with preliminary results from FNAL |p,pr They are 
plotted in Fig. |^ (the Wuppertal-HLRZ group ex- 
trapolated the form factors at = only, using 
various ansatze, so their results are not displayed). 

For massless leptons, the decay rate is determined 
by /+ alone. However, the constraint /^(O) = /"(O) 
(with suitable conventions) makes lattice measure- 
ments of both form factors useful. One procedure 
uses dispersive constraints to obtain model indepen- 
dent information, starting from the two-point corre- 
lation function of = u'^^b, 



W{q^) 



*a;e*9-^ (0| TV^''(a;)t/'^t(o) |o) 



[q^-q^-q^gn^Tiq^] 



Because the contributions of intermediate B* , Bit, 
. . . states to the corresponding spectral functions are 
all positive, the following inequality holds: 

lmUL{t)>m\f(t)\', 

where t = q^ and is a known function. A simi- 
lar inequality relates ImHy and Combining the 
inequalities with subtracted dispersion relations for 
TIt,l provides bounds on the form factors in terms 
of quantities which can be evaluated in perturbative 
QCD. 

Known values of the form factors can be incorpo- 
rated to tighten the bounds. Lellouch [|l^ has shown 
how to incorporate the f~^{0) = /°(0) constraint to- 
gether with imperfectly known values of the form fac- 
tors, typical of lattice results with errors, to obtain 
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Table 2 

Lattices used for calculations of heavy-to-light matrix elements. All simulations are quenched, 
case heavy quarks are treated with the Fermilab formalism 
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families of bounds with varying confidence levels. A 
set of such bounds are shown in Fig. ^ together with 
the UKQCD ||l^ lattice measurements used to obtain 
them. In the figure, /" and /"*" are plotted back-to- 
back, showing the effect of imposing the constraint 
at §2 = 0. 

In Table || the bounds have been used to give 
ranges of values for the decay rate T{B^^tt^I~Di) 
in units of |14bpps~^ together with values for the 
form factor /+ at = 0. When combined with the 
experimental result for the decay rate in Eq. (|l|) , one 
can extract \Vub\ with about 35% theoretical error. 
Although this result is not very precise, the proce- 
dure used relies only on lattice calculations of ma- 
trix elements and heavy quark symmetry, together 
with perturbative QCD and analyticity properties 
in applying the dispersive constraints. There is no 
model dependence. Development and application of 
improved lattice actions for heavy quarks may in fu- 
ture remove the need for heavy quark symmetry in 
the extrapolation from D to B mesons. It is tempt- 
ing to consider that lattice results could eventually 
supplant models in the experimental efficiency deter- 
minations, thereby removing model dependence from 
the experimental results. 

Also shown in Table ^ for comparison, are values 
obtained from lattice calculations where assumed 
dependences have been imposed. For ELC |Q and 
APE ||l^ , one value of /+ has been used, at the given 
value of , fitted to a single pole form with pole mass 
rup. The UKQCD result |jl^ is obtained from a com- 
bined dipole/pole fit to all measured / points. 
Note that the UKQCD points have statistical errors 
only and have not been chirally extrapolated — they 
correspond to a pion mass of around 800 MeV (a sim- 
ilar caveat applies for the FNAL |^,|| results for /+ 



Table 3 

Results for B^-^n^l^vi from dispersive constraints 
applied to lattice results ||l^, together with results 
obtained using ansatze for the form factor /+. Col- 
laboration labels in the left hand column refer to 
lattices described in Table H. The decay rates are 
values for the quantity T{B -^'iT^l~i>i)/\Vub\'^'psr^ 





Rate 


/+(0) 


Dispersive 
Constraint 




2.4-28 
3.6-17 
4.8-10 


-0.26-0.92 95% CL 
0.00-0.68 70% CL 
0.18-0.49 30% CL 


ELC 9±6 0.10-0.49 
q2~18GeV^ pole fit, TOp=5.29(l) GeV 


APE 8 ± 4 0.23-0.43 
(J2~20.4GeV^ pole fit, mp=5. 32(1) GeV (l| 


UKQCDa 7 ± 1 
dipole/pole fit to /+ 


0.21-0.27 



and J^). The results given have used these values 
as though they applied to the physical pion. In ob- 
taining bounds based on these points, Lellouch |]l^ 
added a conservatively estimated systematic error in- 
cluding terms to account for the chiral extrapolation 
(this error has been added to the UKQCD and FNAL 
points plotted in Fig. ||). Needless to say, improved 
lattice results can be used as input for the bounds 
once they become available. 

Discretisation errors are an important source of 
systematic uncertainty Mass dependent er- 

rors affect the crucial extrapolation from D to B 
mesons in calculations using standard Wilson or 
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Figure 2. Lattice measurements of the form factors for i?°-^7r+Z vi and B'^—^p^l vi. Points are from ELC [Q 
(ee), ape (=), UKQCD (•) and FNAL §,| (x, preliminary). Systematic errors have been added 

to the UKQCD and FNAL points, following |n^. The vertical dashed lines show 
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Sheikholeslami-Wohlert (SW) quarks. Momentum 
dependent discretisation errors and increasing statis- 
tical uncertainty also limit the range of spatial mo- 
menta and hence q^ that can be probed. 

Discretisation errors in the heavy-mass scaling 
were studied by the FNAL group using the 

Fermilab improvement formalism. They studied 
(ttI \P) /\/2ME^, where P is a pseudoscalar me- 
son of mass M . This quantity should be constant in 
the M ^ oo limit. The matrix elements were com- 
puted for D and B mesons as well as in the static 
limit. The improvement prescription resulted in lit- 
tle dependence on 1 /M for a range of pion momenta. 

Momentum-dependent discretisation errors 
have been studied by a FNAL-IUinois-Hiroshima 



group |20[| also using the Fermilab formalism for 
heavy quarks. From a study of the axial current ma- 
trix element between a pseudoscalar and the vacuum, 
they estimate that momentum dependent errors are 
less than 20% for light mesons with |p| < 1.2 GeV 
on a relatively coarse P — 5.7 lattice. 

Improvement will clearly help in the determination 
of reliable phenomenological results from lattice cal- 
culations. It will be particularly interesting to test 
the benefit of the 0{a) improvement program out- 
lined by the ALPHA coUaboration 

I note finally in this section that chiral extrapola- 
tions are severe for the _B — s- tt matrix elements. As 
the pion mass approaches its physical value, the B*- 
pole and the beginning of the Bn continuum move 
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GeV^ 

Figure 3. Bounds on the form factors and 
/° for B'^ —>Tr~^ l~ Di obtained from dispersive con- 
straints |17|. The points displayed are from 
UKQCD |1C], with added systematic errors, and the 
shaded band is the prediction of a Hght-cone sum 
rule calculation [03. 
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Figure 4. UKQCD calculation of the differential de- 
cay rate for B^—^p^l^Pi as a function of 



very near the upper endpoint of q^. The form factors 
may vary rapidly with the pion mass near ^i^ax- For 
this reason, the lattice calculations of semileptonic 
B —f p decay are currently most reliable and I now 
turn to this. 

2.2. Semileptonic B ^ p 

To avoid models for the q^ dependence of the form 
factors, we use the lattice results directly. The lat- 
tice can give the differential decay rate dV/dq^, or 
the partially integrated rate, in a region near 
g^ax up to the unknown factor |Vuf,p. For exam- 
ple, UKQCD parametrised the differential decay 
rate near g^^x by, 

^ = const|Kfapg^AV2a2(l -^ (6) 

where A is the usual phase space factor and a and h 
are constants. The constant a plays the role of the 
Isgur-Wise function evaluated at w = 1 for heavy-to- 
heavy transitions, but in this case there is no sym- 
metry to determine its value at leading order in the 
heavy quark effective theory. 

For massless leptons, the differential decay rate de- 
pends on the V , Ai and A2 form factors, but Ai is 
the dominant contribution near g^^x ^^'^ ^^^o the 
best measured on the lattice, as shown in Fig. ^. 

Fig. ^ shows UKQCD results for the differential 
decay rate together with a fit to the parametrisa- 
tion in Eq. @. Because of phase space suppression, 
the point at g^^x does not influence the fit, and so 
gives another determination of a^. The values found 



are 0, 




21(3) GeV^ fit to 5 points, 
23(2) GeV' gLx Point, 



leading to: 

a = 4.6l:[5;|(stat) ± 0.6(syst) GeV. 

Discounting experimental errors, this result will al- 
low determination of \Vub\ with a theoretical uncer- 
tainty of 10% statistical and 12% systematic. With 
full reconstruction of events, CLEO are beginning to 
extract the differential decay distributions The 
Babar experiment, with unequal beam energies al- 
lowing verification that the p and lepton originate 
from the same vertex, should do even better. 

The UKQCD results for V, Ai and A2 agree very 
well with a light cone sum rule (LCSR) calculation 
of Ball and Braun ||2|]. The agreement is perhaps 
fortuitous. More interestingly, LCSR calculations 
predict that all the form factors for heavy-to-light 
decays have the following heavy mass dependence at 
q^ = (I9|j2|]: 

/(O) = + + + •••)• (7) 

The leading M dependence comprises \/M from 
the heavy state normalisation together with the be- 
haviour of the leading twist light cone wavefunction. 
For Ai{0) the LCSR result is || Ai{0) ~ 0.26. 
UKQCD fitted the heavy mass dependence of ^i(O) 
to the form in Eq. |^ and found ijlll 

r 0.18 ± 0.02 ao and ai ^ 
^i^'^^ - \ 0.22l[;;!]4 ao, ai and ^ 
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Table 4 

Lattice results for B^K*j. The labels in the left 
hand column refer to the lattices described in Ta- 
bleg. 
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UKQCDa 
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0.30(3) 



Pole fits for Ai (and dipole/pole for / f'^ for _B tt 
decays) have leading M~'^/^ behaviour at q-^ = 0, so 
we can also compare with other lattice results: 

r 0.22 ±0.05 ELC ]ll 
Ai(0) = <^ 0.24 ±0.12 APE|5| 

[ 0.27tl UKQCD 

The Wuppertal-HLRZ group |l^,|l^ tried various 
forms for the heavy mass dependence of Ai{0), al- 
though none had a leading Af^"^/^ dependence. 



2.3. Rare radiative B K*j 

This decay was discussed in some detail by A. Soni 
at Lattice 95 |Q so my comments will be brief. In 
Table ^ I summarise the available lattice results, all 
from quenched simulations for the matrix element 



(8) 



which is parameterised by three form factors, T^, 
i = 1,2,3. For the decay rate the related values 
Ti(0) and 12(0) are needed. Suitably defined, they 
are equal, so in the table I quote a single value T(0), 
together with the directly measured T2{q^^y-). The 
results are classified according to the leading M de- 
pendence of the form factor at = which is gov- 
erned by the model used to fit the q^ dependence. 
Dipole/polc forms for T1/T2 give M~^/^ behaviour 
and pole/constant forms give The table 

shows that the results agree when the same assump- 
tions are made. All groups find that T2 has much 
less dependence than Ti, but the overall forms 
cannot be decided, so a phcnomcnological prediction 
is elusive. 

Additional long distance contributions may not be 
negligible so the matrix element of Eq. (||) may not 
give the true decay rate |p5|-|27[| . Once the depen- 
dence of the form factors is known, lattice calcula- 
tions of the ratio Rk' = T{B^K*j)/r{b 57) can 
be compared to the experimental result in Eq. (H) to 
test for long distance effects. 



Heavy quark symmetry, combined with light 
flavour SU (3) symmetry, relates the form factors 
for B°^p+l-i?i and B^K*-f in the infinite 

heavy quark mass limit. On the lattice these re- 
lations can be tested using identical light quarks. 
UKQCD showed that the ratios V/2Ti and 
A1/2T2 both satisfied the heavy quark symmetry 
prediction of unity in the infinite mass limit. A com- 
bined fit of the pseudoscalar to vector form factors 
consistent with heavy quark symmetry could help 
resolve the ambiguity in the q^ dependence of the 
B^K*^ form factors. 

3. LEPTONIC DECAY CONSTANTS OF 
HEAVY-LIGHT PSEUDOSCALARS 

Leptonic decay constants of heavy light systems 
were fully reviewed by C. Allton at Lattice 95 |30|, so 
here I report new results. To establish notation, the 
decay constant is determined by lattice calculation 
of the dimensionless quantity Zl according to, 

where Z'''" is the renormalisation constant required 
to match to the continuum and a is the lattice spac- 
ing. 

3.1. Conventional Methods 

New results were shown by the MILC and 
JLQCD |3^ collaborations. Both groups perform 
continuum extrapolations from results at several /3's. 

MILC have six sets of quenched configurations for 
(3 in the range 5.7-6.52, together with six sets of two- 
flavour dynamical staggered fermion configurations 
with 5.445 < /3 < 5.7. A Wilson action is used for 
the valence quarks, with a hopping parameter expan- 
sion [ p3[ for the heavy quark for easy simulation of 
a range of masses. The heavy quark has a Kronf eld- 
Lepage-Mackenzie normalisation, \/l ~ 6k, together 
with a shift from the pole mass to the kinetic mass 
for the heavy meson. The continuum extrapolation 
should deal with remaining 0{a) errors. Fig. ^ shows 
this extrapolation for Jb- MlLC's (preliminary) re- 



suits are: 




fB = 


166(11)(28)(14) MeV 


/b. = 


181(10)(36)(18) MeV 


IBs I fB = 


1.10(2)(5)(7) 


Id = 


196(9)(14)(8)MeV 


/d. = 


211(7)(25)(ll)MeV 


Id, I fo = 


1.09(2)(5)(5) 



The central values come from the quenched results 
using a linear extrapolation in a with the scale set 
by /tt and linear chiral extrapolations (in 1/k). The 
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Figure 5. Continuum extrapolation of Jb by 
MILC The dashed line is a constant fit for 

quenched results with /? > 6.0. 

heavy mass dependence is determined by a fit for 
heavy meson masses in the range 1.5 GeV < M < 
4 GeV combined with a point in the static limit. The 
first error quoted is a combination of the statistical 
error and that coming from choice of fitting proce- 
dure, the second error is remaining systematic errors 
within the quenched approximation and the final er- 
ror is for quenching. 

Improved covariant fitting methods have changed 
the results compared to their values at Lattice 
95 [3^ ], although within errors. In particular the 
difference between quenched and unquenched results 
is less dramatic, but still looks significant (especially 
for Jb^i and indicates an increase in the value of the 
decay constant for full QCD. 

JLQCD have results from three /3 values, 5.9, 
6.1 and 6.3 with quenched configurations using the 
Wilson action for both light and heavy quarks. They 
study different prescriptions for reducing the 0{ma) 
scaling violations associated with the heavy quark 
and aim to show that results from all prescriptions 
converge in the continuum limit. Smeared wavefunc- 
tions are determined for each available combination 
of heavy and light quark kappa values. Results are 
quoted using the charmonium IS- IP splitting to set 
the scale. 

JLQCD apply four different procedures to extract 
the decay constant. The first uses the traditional 
\^2k quark field normalisation and the meson pole 
mass. Three further methods use a Kronfeld-Lepage- 
Mackenzie normalisation combined with, (a) setting 
the meson mass from the pole mass, (b) pole mass 
shifted to kinetic mass at tree level, or (c) measured 
kinetic mass. For all four methods, a linear extrapo- 



Figure 6. Continuum extrapolation of /s by 
JLQCD using different prescriptions. 

lation in a gives the continuum result. This is shown 
for Jb in Fig. ^. The preliminary results are js^ : 

fB = 163(19)(i«)MeV 
fB^ = 174(12) (i^)MeV 

fo = 191(8)(^)MeV 
Jd, = 210(6) (il)MeV 

where the first error is statistical and the second is 
from the spread over the four methods above. 

The MILC and JLQCD results are included in 
Table | which summarises determinations of heavy- 
light pseudoscalar meson decay constants and their 
ratios. MILC and JLQCD are (currently) in pleas- 
ing agreement. From the table, and allowing for an 
increase in the value from unquenching, I quote a 
global result for Jb- 

fB = 175± 25MeV. 

3.2. Decay Constants from NRQCD 

The SCO collaboration have calculated matrix el- 
ements for determining pseudoscalar and vector B 
meson decay constants [^,^ using both Wilson and 
tadpole-improved SW light quarks combined with 
NRQCD heavy quarks. They have 100 unquenched 
configurations with nj — 2 flavours of dynamical 
staggered quark {mscaO, — 0.01) at /3 = 5.6. The 
NRQCD action is corrected to 0{1/M) and matrix 
elements are calculated for extra operators appear- 
ing at 0{1/M) in the matching of NRQCD currents 
onto full QCD vector and axial vector currents^. 

Appropriate renormalisation constants are not yet 
available, however, for the matching to full QCD, so I 

^See Refs. |45y46[ for details and qualifications. 
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Table 5 

Summary of results for values and ratios of leptonic decay constants of pseudoscalar mesons. Statistical and 
systematic errors have been combined in quadrature. The results from MILC | |3l[ | include a systematic error for 
quenching. Results from APE and FNAL show statistical errors only. FNAL ||] results are preliminary 
and use a nominal value for the strange quark mass to determine Jb, and fo^ ■ Results with no quoted values 
for (3 or come from more than one lattice: PSI-WUP [Q, MILC |Q and JLQCD perform a continuum 
extrapolation using their own data. 

Ref (3 



FNAL 

DeG-L 

APE 

UKQCD 

LANL 

BDHS 

UKQCD 

BLS 

ELC 

ELC 

ELC 

PSLWUP 
LANL 

PSLWUP 

MILC 

JLQCD 



111 

137] 
135] 
|3S1 
139] 

m 

138] 
141] 

11 

ll 

142] 

m 
m 

m 

131] 

ll 



/3 




fB 






Id 






GeV 


MeV 


MeV 


Ib 


MeV 


MeV 


Id 


5.9 




188(1) 


207(1) 




220(1) 


239(1) 




6.0 


1.9 








190(33) 


222(16) 


1.17(22) 


6.0 




197(18) 






218(9) 


240(9) 


1.11(1) 


6.0 


2.0 


mil) 




1.17(12) 


199(11) 




1-13(?) 


6.0 


2.33 








229(?i) 


mil) 


1.14(2) 


6.1 


2 








174(53) 


234(72) 


1.35(22) 


6.2 


2.7 






1.22(1) 


185(1) 


212(1) 


1.18(2) 


6.3 


3.0 


187(38) 


207(41) 


1.11(6) 


208(38) 


230(36) 


1.11(6) 


6.4 


3.3 








210(40) 


230(50) 





205(40) 1.08(6) 



180(50) 1.09(5) 
166(33) 181(41) 1.10(9) 
mil) 174(i) 



194(15) 

198(17) 209(18) 

186(29) 218(15) 

170(30) 

196(18) 211(28) 

191(i?) 210(ii) 



1.09(5) 
1.09(7) 



do not quote here any values for decay constants. Ac- 
cepting this limitation, the slope of the pseudoscalar 
decay constant with respect to 1/M is calculated, 
that is, the term cp in the expansion 

fy/M = const(l + cp/M + •••)■ 

The results are, 




-1.35(15)a-i = -2.8(5) GeV Wilson, 
-1.0(2)a-i - -2 GeV TI SW, 



to be compared to slopes of around 1 GeV found us- 
ing conventional methods. Study of a spin average 
of pseudoscalar and (appropriately defined) vector 
meson decay constants shows that just over 90% of 
the effect comes from the heavy quark kinetic en- 
ergy term. Analysing the effect of the kinetic term on 
wavefunctions and energies in the static limit (from a 
spinless relativistic quark model) to first order shows 
that a slope of 1-2 GeV is natural, and also predicts 
that the IS state energy should rise slowly with 1 /M, 
with the 2S energy rising faster 

It will be very interesting to see the slope and de- 
cay constant values once the correct renormalisation 
constants are included. 



3.3. f§^* from Bermions 

The APETOV group have results for /|*'^* ob- 
tained using a pseudofermion method [4^ ]. In stan- 
dard static quark methods, smeared sources must 
be tuned to optimise the projection onto the ground 
state and avoid noise. Moreover, the numerical ex- 
pense of inverting the Dirac operator restricts the 
set of available light quark propagators to have one 
endpoint fixed. In contrast, APETOV determine the 
light propagator by Monte Carlo inversion, allowing 
them to average the static correlator over all points 
without requiring smeared sources. The additional 
average over lattice points helps overcome the extra 
noise inherent in the Monte Carlo inversion. Further 
averaging over the gauge links is also employed. 

The bermion action is J2x \Q4>{x)\'^, where x is 
a lattice point, U is the gauge field, (j) the scalar 
bermion field and Q is given by Q4>{x) — j5D(l>{x), 
where D is the usual Wilson lattice Dirac operator. 
The bermions are thermalised in each gauge config- 
uration and the pseudofermion two-point function 
then gives (Q^)"^. Rcmultiplying by Q and 75 al- 
lows the usual light quark Dirac propagator to be 
determined between any lattice points. The corre- 
lator of two local heavy-light bilinears with a given 
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time separation is then averaged over all points to 
determine Zj^. In practice, a further correlator using 
a (Q'^)"^ bermion propagator with the static prop- 
agator, which better isolates the lowest lying state, 
needs to be analysed concurrently. 

The method has been applied to two 16'^ x 32 lat- 
tices with 30 gauge configurations each at /3 = 5.7 
and 6.0, with the results B^: 



0.477(31) Kent 

0.582(11) 

^strange 
0.188(37) Kerit 

0.216(13) 

^stranec 



^ = 5.7 



(i = 6.0 



The /? = 6.0 chirally extrapolated result agrees well 
with world results for Z^. At /3 = 5.7 the result is 
below the FNAL Q value Zl = 0.564(28), but is 
on a larger lattice and has comparable errors with 
one third the number of gauge configurations. 

APETOV have extended the method by making 
the bermions dynamical |51|. Each bermion flavour 
corresponds to —2 fermion flavours. Zj, is calculated 
for different n;, = —2nf, tuning /3 to keep the ratio 
R = Tn^/rn^p constant, starting with fi = 5.7 for 
rib = 0. Including a one-loop perturbative value for 
Z''®'^, the result for the ratio of Z^Z''"" for 3 flavours 
to flavors is 1.14(2) for R = 0.6 and 1.16(4) for 
R = 0.5, with the effect increasing to about 20% in 
the chiral limit. This offers more evidence for the 
increase of Jb once dynamical effects are included. 

4. Results for B Mixing 

Recent results for the B meson mixing parame- 
ter, Bb(m) _= {B\0{y) \B) /{8/3)flMl, where O = 
67^(1— 75)967^(1— 75)9 is the AB = 2 operator and 
Mb is the B meson mass come from UKQCD [Q, 
APE |5^ and the Kentucky group |Q in the static 
limit, together with results ||2^,^ using conventional 
methods. The APE result makes use of a new cal- 
culation |5^ of the full-theory/static-theory match- 
ing which incorporates previously omitted contribu- 
tions. Table ^ is a compilation of lattice determina- 
tions of Bb- To ease comparison, values are given for 
-B_b('ti{,=5 GeV) and for the (1-loop) renormalisation 
group invariant quantity, Bb = o:s{i^)^'^^^° Bb{i^)- 
A similar collection of results appears in the contri- 
bution by Christensen, Draper and McNeile |^ to 
these proceedings. 

The conventional results are consistent and show 
no a dependence. The static results show substan- 
tial differences arising from the renormalisation con- 
stants. The Bb calculation involves dividing by the 
square of Za, the heavy- light axial current renor- 
malisation, which differs for Wilson (KEN) and SW 
(UKQCD, APE) light quarks. Moreover, whether 



products of renormalisation constants are expanded 
to a given order in a or simply multiplied makes a 
significant difference: the UKQCD result for Bb rises 
to 1.19(6) [|4j when the factors are multiplied rather 
than expanded. Nonperturbative determinations of 
the renormalisation factors are clearly crucial to re- 
duce systematic errors ]59|j60[ . 

The relevant phenomenological quantity is /b^b 
which can be extracted directly from the AB = 2 
matrix element. To avoid uncertainties from set- 
ting the scale, it is convenient to determine the ra- 
tio fg^BBjf%BB- BbJBb is found to be close 
to unity in recent calculations: 1.01(3)(3) pil and 
1.011(8) ||. For /s,//b, the values in Tableg can 
be combined with static results. 



f stat 

f stat 
J B 



1.11(2) 
l-13(t) 







= <; 1.22(4)(2) [50| 
[54| 
[53 



1-16(^) 
1.17(3) 



61| 



The second error in the first result iox BbJ Bb is for 
quenching based on numerical evidence for a small 
increase in the ratio on nf = 2 dynamical configura- 
tions (Sharpe and Zhang | |6^ , |63t estimate a quench- 
ing error of —0.04 based on chiral loops). Unquench- 
ing is expected to increase the value of fB^/fs by 
about 10% (the chiral loop estimate in 0.16 [ p2p^ ]). 

Bernard, Blum and Soni | |55| , |6^ report a prelimi- 
nary value for 



Ib^b 



1.81(8)(25). 



^BJB 



The ratio shows no evidence of a dependence on a 
range of lattice spacings. The result implies rather 
a large value for /b^/Zb, compared to other lattice 
results: it will be interesting to see the value of this 
ratio when the decay constants are extracted sepa- 
rately from the same data. 

5. OTHER RESULTS 

Heavy-to-heavy quark transitions have been omit- 
ted from this report. See j6|] for a recent com- 
pilation of lattice results for the Isgur-Wise func- 
tion relevant for B D^*'> semileptonic decays. 
UKQCD presented preliminary results at this confer- 
ence for the baryonic Isgur-Wise function in semilep- 
tonic A;, — + Ac decays Q. 

Results for A2 , the matrix element of the chromo- 
magnetic moment operator between heavy mesons 
in the heavy quark effective theory, continue to sit 
at about half the experimental value The 
one loop renormalisation constant is uncomfortably 
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Table 6 

Lattice results for the meson mixing parameter Bb- Lattice parameters are quoted except where the 
authors have performed a continuum extrapolation from several results, indicated by a— s-0 in the column. 
The column labelled "Action" shows the action used for the light and heavy quarks respectively. The b quark 
mass TOfc has been set to 5 GeV: where a result has been quoted at a different scale, the scale fi and BsifJ-) 
have been listed. The results headed "Static fitted" are obtained by extrapolating 1/M — > where M is the 
heavy-light meson mass. Results in oblique type have been obtained from the authors' values using one-loop 
scaling with 5 flavours and A^.^^^ — 130 MeV. 



Rcf 



/3 



Cfgs 



GeV 



Action 



GeV 



;(^) BB{mb) 



Static 

KEN H 

APE 153] 

UKQCD II 



6.0 20^ X 30 
6.0 243 X 40 
6.2 243 X 48 



32 
600 
60 



2.1 
2.0 
2.9 



W-stat 

SW-stat 

SW-stat 



4.33 0.98(5) 



0.97(5) 
0.82(4) 
0.69(4) 



1.45(8) 
1.21(6) 
1.02(6) 



Static fitted 
ELG jll 

ig 



BS 



6.4 243 X 60 
5.7-6.3 



20 



3.7 

a^O 



W-W 
W-W 



3.7 

2 



0.90(5) 
1.04(5) 



0.88(5) 
0.96(5) 



1.30(7) 
1.42(7) 



Conventional 



BBS [ 


55 


,56 


1 5.7 


163 


X 


33 


100 


1.45 


W- 


-w 


2 


0.96(5) 


0.89(5) 


1.32(7) 


BBS 


55 


,56 


1 6.0 


163 


X 


39 


60 


2.06 


W- 


-w 


2 


0.96(5) 


0.89(5) 


1.32(7) 


BBS 


55 


,56 


1 6.0 


243 


X 


39 


40 


2.22 


w- 


-w 


2 


0.98(4) 


0.90(3) 


1.34(5) 


JLQCD 


57 


] 6.1 


243 


X 


64 


200 


2.56 


w- 


-w 






0.90(5) 


1.32(7) 


JLQCD 


57 


6.3 


323 


X 


80 


100 


3.38 


w~ 


-w 






0.84(6) 


1.24(9) 


BBS 


55 


@ 6.3 


243 


X 


61 


60 


3.40 


w- 


-w 


2 


1.09(16) 


1.01(15) 


1.49(22) 


ELG 


42 


6.4 


243 


X 


60 


20 


3.7 


w- 


-w 


3.7 


0.86(5) 


0.84(5) 


1.24(7) 


BDHS 


40 


] 5.7- 


-6.1 








a^O 


w~ 


-w 


2 


1.10(15) 


0.94(14) 


1.38(21) 


BS 


24 


] 5.7- 


-6.3 








a^O 


w- 


-w 


2 


0.96(6) 


0.89(6) 


1.31(8) 



large: this and other systematic effects warrant fur- 
ther study given the importance of A2 in inclusive B 
decays. 

Semileptonic D — > if *^*) decays were reviewed by 
J. Simone at Lattice 95 1^. 
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